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Abstract 


The  U.S.  Navy  is  interested  in  investigating  the  use  of  ocean  color  data  for  naval  applications,  as  they 
apply  within  the  Air  Defense  Initiative.  The  major  source  of  available  data  for  assessment  is  the  Coastal  Zone 
Color  Scanner  (CZCS)  data  base.  These  data  were  collected  during  1979-1986.  Although  the  CZCS  is  no 
longer  operational,  much  information  has  been  published  on  the  analysis  of  the  data.  This  technical  note 
reviews  that  literature,  which  deals  primarily  with  the  measured  and  estimated  errors  associated  with  using 
this  CZCS  data  to  derive  either  chlorophyll  concentrations  (C)  or  the  diffuse  attenuation  coefficient  (K)  of 
oceanic  waters. 

The  recognized  uncertainties  in  the  absolute  values  of  the  quantities  of  C  and  K  must  be  taken  into 
account  when  addressing  the  potential  use  of  CZCS  data  for  certain  naval  applications.  The  literature 
generally  agrees  that  on  a  global  basis,  K  or  C  can  be  determined  from  CZCS  data  to  only  within  a  factor 
of  2,  as  compared  with  in  situ  measurements,  in  addition,  the  water-leaving  radiance  for  one  spectral  band 
can  probably  be  obtained  from  the  CZCS  data,  but  in  a  more  regional  or  smaller  geographical  area,  after 
corrections  to  within  +.15%.  This  percentage  relates  to  an  error  of  approximately  +30%  in  K  or  C,  both  of 
which  use  a  ratio  of  at  least  two  spectral  bands.  Extreme  care  must  be  taken  in  applying  the  appropriate 
corrections  on  a  pixel-to-pixel  basis  if  the  error  is  to  be  reduced  to  this  value.  If  the  concern  is  with  only 
spatial  trends  of  water  clarity  at  a  particular  point  in  time  or  with  temporal  trends  at  each  point  in  space,  the 
CZCS  images  can  probably  be  used;  however,  the  conditions  under  which  they  apply  must  be  indicated. 

A  preliminary  computer  simulation  that  focuses  on  the  impact  of  some  of  the  data  errors  is  also  included 
with  this  review.  Further  modeling  could  lead  to  improved  atmospheric  corrections,  leading  to  required 
statistical  variances  of  the  pixel  radiance  values  as  a  function  of  geographical  region  and  spectral  band. 
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I.  Introduction 


The  Coastal  Zone  Color  Scanner  (CZCS),  flown  on  the  NIMBUS-7  satellite,  collected  ocean  color  data 
from  July  1979  to  June  1986.  This  scanning  radiometer  viewed  the  ocean  in  six  coregistered  spectral  bands; 
five  in  the  reflective  region  of  the  spectrum  (centered  at  443,  520,  550,  670,  and  750  nm)  and  one  in  the 
thermal  IR  (10.5  to  12.5  jim).  Although  the  CZCS  data  has  contributed  enormously  to  our  knowledge  of  the 
world’s  oceanic  productivity,  its  data  is  known  to  include  errors.  Many  of  these  errors  are  associated  with 
radiometric  sensitivity  of  the  instrument,  exact  derivation  of  the  pigment  concentration,  the  effect  of  subpixei 
sized  clouds  on  the  retrieved  radiances,  complete  removal  of  the  atmospheric  turbidity,  and  navigational 
errors  to  cite  a  few. 

This  technical  note  summarizes  the  results  of  a  number  of  key  technical  reports  concerned  with  the 
errors  associated  with  comparing  in  situ  ocean  color  measurements  with  those  derived  from  satellite  CZCS 
data.  A  fairly  comprehensive  bibliography  has  been  included. 

A  preliminary  computer  based  sensor  error  model  was  constructed  and  presented  in  section  IV.  It  is 
given  to  illustrate  the  potential  value  of  using  a  modeling  approach  to  assess  the  nature  and  impact  of 
various  environmental  parameters  on  the  total  error  budget  of  the  CZCS  data. 

II.  Satellite  Derived  Ocean  Color  Data 

A.  Algorithms  for  correlating  satellite  and  in  situ  measurements 

“Ocean  Color"  is  a  measure  of  the  spectral  reflectance  of  the  water  column.  It  is  measured  remotely  as 
the  wavelength-dependent  radiance  emerging  from  the  sea  surface.  Ratios  of  water  leaving  radiance  are 
used  to  calculate  phytoplankton  pigment  concentration.  Research  using  CZCS  data  has  established  the 
feasibility  of  using  ocean  color  observations  from  space  to  study  global  and  mesoscale  distributions  of 
ocean  bio-optical  properties.  The  CZCS  data  has  been  used  extensively  for  its  potential  ability  to  map  global 
distributions  and  standing  stock  of  marine  phytoplankton. 

The  "diffuse  attenuation  coefficient"  (K(X))is  a  quantity  that  has  been  extensively  used  as  a  measure  of 
water  quality  or  visibility.  This  quantity  is  of  primary  importance  to  those  naval  applications  that  are 
concerned  with  operation  of  active  and  passive  optical  sensors. 

Smith  and  Baker  (1978;  1982)  related  the  spectral  diffuse  attenuation  coefficient  (K(X))  to  chlorophyll 
pigment  concentrations  (C),  while  Clark  (1981)  derived  an  algorithm  relating  (CX),  to  the  upwelling  radiances 
Lw(443)  and  Lw(550)  given  below 
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Austin  and  Petzold  (1981)  in  their  classic  article  on  this  subject  related  the  upwelling  radiances  at  Lw(443) 
and  Lw(550),  and  (C)  through  the  following  empirical  expressions 
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where: 

Lw(443)  =  inherent  upwelling  radiance  at  443  nm  at  the  ocean  surface 
Lw(550)  =  inherent  upwelling  radiance  at  550  nm  at  the  ocean  surface 
CX  =  average  chlorophyll  pigment  concentration  in  mg  mV3. 

These  algorithms  were  developed  using  CZCS  data  and  in  situ  spectroradiometric  data  from  88 
oceanographic  stations  from  a  large  number  of  investigators  working  In  a  wide  variety  of  water  types.  They 
concluded  that  it  should  be  possible,  under  most  field  conditions,  to  infer  values  of  C  or  K  from  satellite 
ocean  color  measurements  with  accuracies  comparable  with  that  of  in  situ  measurements.  They  also  stated 
that  this  was  very  dependent  upon  one’s  ability  to  make  the  appropriate  atmospheric  corrections. 

B.  Correlation  of  satellite  derived  and  in  situ  surface  chlorophyll  concentrations 

Gordon  et  al.  (1980)  compared  CZCS  and  surface  distributions  and  in  situ  pigment  concentrations  in 
the  Gulf  of  Mexico  during  a  Nimbus-7  post  launch  experiment  in  November  1979.  They  found  that  in  the 

worst  case  the  pigment  concentration  (C)  could  be  estimated  from  radiance  ratios  to  within  a  factor  of  2, 

and  the  ratio  of  the  radiances  of  Lw(443)/Lw(550)  was  useful  only  for  pigment  concentrations  of  less  than 
0.6  mg  m-3  (K= 0.089-equation  3).  A  number  of  papers  appeared  as  a  result  of  the  South  African  Ocean 
Color  and  Upwelling  Experiment  for  the  time  period  November  1978  through  April  1980.  (Shannon  et  al., 
1985;  Walters,  1985;  and  Walters  et  al.,  1985)  The  results  of  these  studies  showed  that  chlorophyll 
concentrations  can  be  derived  from  CZCS  radiance  data  with  an  accuracy  of  better  than  a  factor  of  2. 

Muller-Karger  et  al.  (1990)  examined  21  CZCS  images  of  the  Southeastern  Bering  Sea  to  map  the  near- 
surface  distribution  of  phytoplankton  during  1979  and  compared  the  results  to  in  situ  measurements  made 
from  the  ship.  They  found  that  the  CZCS  data  underestimated  the  ship  data  in  chlorophyll  pigment 
concentration  by  a  factor  of  at  least  2  during  the  spring  of  1979.  During  a  similar  period  of  time  in  1980  the 
two  sets  of  data  agreed  much  better.  They  concluded  that  region  specific  pigment  algorithms  may  be 
required  to  yield  quantitative  results  from  remote  sensing  data — at  least  from  the  Bering  Sea.  Similarly,  even 
at  more  moderate  latitudes,  the  sun  elevation  may  be  sufficiently  low  In  the  winter  to  give  rise  to  large  errors 
in  the  inferred  values  of  pigment  concentrations. 

Gordon  et  al.  (1983a)  compared  CZCS  data  with  ship  data  taken  aboard  the  RV  Athena  II NOAA  -  NESS 
CZCS  ocean  color  cruise  (May  31,  1979  -  June  23,  1979)  from  the  Gulf  of  Maine  to  Florida.  The  results 
showed  that  the  CZCS  derived  values  for  C  reflected  a  ±30-40%  difference  with  in  situ  values  of  C  within  the 
pigment  concentration  range  0.08  -1.5  mg  m-3.  These  results  were  obtained  for  regions  showing  relatively 
weak  horizontal  gradients,  i.e.,  horizontal  scale  of  variability  significantly  greater  than  the  CZCS  pixel  size. 
A  total  of  four  images  of  the  shelf  and  slope  waters  of  the  Middle  Atlantic  Bight  and  the  Sargasso  Sea  were 
processed  and  averaged.  The  satellite  derived  values  of  C  were  compared  with  in  situ  values  obtained  from 
continuous  measurements  made  from  ships.  The  quoted  errors  are  rms  values. 

Gordon  and  Clark  (1981)  presented  data  comparing  in  situ  pigment  concentrations  and  upwelling 
radiance  for  60  Case  I  water  sites  off  the  west  and  east  coasts  of  the  U.S.  and  the  Gulf  of  Mexico.  (Case 
I  waters  are  defined  as  those  waters  for  which  phytoplankton  and  their  detrital  material  play  the  dominant 
role  in  determining  the  optical  properties.  The  majority  of  open  ocean  waters  can  be  classed  as  Case  I.  The 
results  showed  that  a  ±20%  error  existed  in  the  upwelling  radiance  (In  one  band)  which  results  in  an 
approximate  error  In  the  pigment  concentration  (C)  of  ±40%  for  concentrations  less  than  0.25  mg  m3. 
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C.  Determination  of  water  clarity  to  depth  via  CZCS  data 

The  above  discussion  has  dealt  mainly  with  the  measurement  of  phytoplankton  pigment  concentrations 
in  the  near  surface  water  column;  however,  in  many  naval  applications  one  is  more  concerned  with  the 
interpretation  of  satellite  data  to  yield  light  penetration  to  depth.  As  shown  in  section  II,  there  is  a  relationship 
between  the  pigment  concentration  (C)  and  water  clarity  (diffuse  attenuation  coefficient).  The  errors 
associated  with  CZCS  derived  optical  data  refer  mainly  to  near-surface  (i.e.,  one  attenuation  length  (1  /K)) 
measurements,  which  in  terms  of  K  values  of  interest  yields  water  depths  ranging  from  10  to  50  m.  However, 
for  many  practical  applications  one  is  concerned  with  the  use  of  these  surface  measurements  to  predict  the 
optical  properties  to  several  attenuation  lengths.  Several  of  the  papers  that  appear  in  the  bibliography  have 
touched  on  this  problem.  Shannon  et  al.  (1985)  sounded  a  warning  about  the  use  of  satellite-derived 
chlorophyll  data  in  the  region  east  of  Cape  Agulhas  during  the  summer  owing  to  the  marked  subsurface 
chlorophyll  maximum  that  is  beneath  the  low  chlorophyll  layer.  They  also  stated  that  one  must  be  careful 
in  extrapolating  CZCS  imagery  to  areas  where  the  algorithms  have  not  been  adequately  validated 

Straus  and  Woods  (1988),  presented  comparisons  of  CZCS  derived  chlorophyll  pigment  concentrations 
with  in  situ  measurements  taken  on  a  series  of  ship  cruises  from  1984  to  1987  between  the  Azores  and 
Greenland.  The  total  data  set  consisted  of  20,000  profiles  (from  10  m  to  200  m).  They  showed  that  only  a 
small  fraction  (5-10%)  of  the  chlorophyll  concentration  is  actually  sensed  by  the  CZCS  in  this  region— mainly 
due  to  the  depth  of  the  chlorophyll  maxima.  They  concluded  that  the  CZCS  derived  chlorophyll 
concentrations  and  those  measured  from  the  ship  differed  by  a  factor  of  2.  There  is  a  requirement  for  both 
experiments  and  theoretical  modeling  to  de'armine  the  impact  of  depth  variation  of  pigment  concentrations 
of  spectral  ieaving  radiance.  The  mear  ,n  situ  chlorophyll  concentration  of  the  euphoric  zone  is 
overestimated  by  a  factor  of  nearly  2  in  early  spring,  when  the  plankton  are  concentrated  near  the  surface, 
and  underestimated  in  late  summer,  when  the  plankton  are  concentrated  in  a  layer  deep  in  the  seasonal 
thermocline.  They  claimed  their  main  source  of  error  (differences  as  large  as  a  factor  of  3  or  4)  arose  from 
the  fact  that  the  algorithms  used  for  the  CZCS  could  not  correct  for  horizontal  and  seasonal  variation  in 
subsurface  chlorophyll — both  in  chlorophyll  content  and  mean  concentration  in  the  euphoric  zone. 

Kitchen  and  Zaneveld  (1990)  examined  relationships  between  beam  attenuation,  absorption,  suspended 
particle  concentrations,  size  distribution,  and  pigment  content.  They  concluded  that  the  remote  sensing 
algorithms  utilized  to  predict  backscattering  from  chlorophyll  concentration  were  developed  using  near¬ 
surface  data  and  should  not  be  used  for  the  prediction  of  vertical  structure  of  backscattering.  Since  the 
CZCS  measures  only  the  upper  one  or  two  attenuation  length  at  best,  layers  below  the  mixed  layer 
contribute  very  little  to  the  CZCS  surface  reflectance.  They  further  stated  that  the  remote  sensing/in  situ 
correlations  should  not  be  used  when  one  is  estimating  the  performance  of  UDAR  systems  or 
bioluminescence  radiance  below  the  mixed  layer.  Other  investigators  such  as  Mueller  and  Lange  (1989) 
indicated  that  the  vertical  optical  structure  in  the  ocean  can  be  predicted  by  the  CZCS  radiance  values  if 
sufficient  prior  knowledge  erf  the  climatology  is  known  along  with  continual  in  situ  sampling  of  the  waters. 
The  1984  California  Cooperative  Oceanic  Fisheries  Investigations  (CalCOFI)  cruise  data  (Thomas  and  Strub 
(1990)  and  Strub  et  al.  (1990))  shows  that  in  general,  the  subsurface  chlorophyll  maximum  north  of  the  front 
is  deeper  than  20  m  over  the  entire  year  and  probably  is  never  seen  by  the  CZCS. 

There  are  numerous  naval  and  oceanographic  applications  that  depend  upon  the  knowledge  of  water 
clarity,  such  as  an  airborne  laser:  for  mapping  underwater  targets,  shallow  water  bathymetry,  alrcraft-to- 
submarine  communications,  (SLCEVAL,  1989a  and  b)  and  for  potentially  measuring  vertical  sound  velocity 
profile  via  Brillouin  scattering  (Hickman  et  al.,  1 991 ).  In  calculating  the  performance  of  such  a  UDAR  system, 
one  of  the  key  environmental  parameters  that  must  be  considered  is  the  water  diffuse  attenuations  coefficient 
(K)  and  its  associated  errors.  As  mentioned  throughout  this  technical  note,  the  best  determination  of  K 
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derived  from  satellite  ocean  color  data  that  one  can  expect  is  approximately  ±30%.  Figure  1  is  given  to 
show  the  effect  that  an  error  of  this  magnitude  in  K  has  upon  the  uncertainty  in  target  depth  detection  for 
an  airborne  UDAR  system  as  a  function  of  the  number  of  attenuation  lengths  and  various  values  of  K.  For 
instance,  assume  a  UOAR  system  has  been  designed  to  have  a  target  detection  depth  (d)  capability  of  4 
attenuation  lengths  (i.e.,  kd =4).  Then,  for  a  nominal  water  clarity  of  k=0.05  m  \  one  would  expect  to  detect 
targets  to  depths  of  80  m.  However,  referring  to  Figure  1 ,  the  uncertainty  of  ±30%  in  K  (at  k=0.05  m'1)  leads 
to  a  depth  uncertainty  of  54  m.  This  results  in  a  value  for  the  predicted  target  depth  capability  for  these 
types  of  waters  of  approximately  80  ±27  m. 


Figure  1.  Total  error  in  depth  of  target  detection  for  an  airborne  laser  system  (The  error  in  K  is 
assumed  to  be  ±30%). 
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III.  Sources  of  Error 


Although  a  great  deal  of  research  has  been  reported  on  the  atmospheric  corrections  that  must  be 
applied  to  the  CZCS  data,  there  has  been  little  mention  of  the  magnitude  of  the  errors  associated  with  other 
corrections  that  should  be  considered.  A  complete  error  analysis  of  the  CZCS  data  has  not  been  made. 

The  errors  associated  with  the  remote  sensing  of  the  optical  clarity  of  the  water  by  the  CZCS,  as  listed 
by  the  various  CZCS  investigators,  arise  from  a  number  of  factors  such  as: 

a.  atmospheric  correction  (single  scattering  or  multiple  scattering) 

b.  sensor  degradation  -  changing  calibration 

c.  subpixel  sized  clouds 

d.  sensor  ringing 

e.  geographical  navigation  error  (ship  and  satellite) 

f.  white  caps 

g.  in  situ  errors  in  determination  of  chlorophyll/pigment  concentration 

h.  contamination  of  the  signal  by  nonbiological  material. 

In  addition  to  the  above  corrections,  polarization  effects  have  been  ignored.  Recent  work  by  Kattawar 
and  Adams  (1989)  indicates  that  it  may  now  be  possible  to  take  polarizations  effects  into  account  and  make 
appropriate  corrections.  The  importance  of  polarization  was  discussed  in  a  recent  report  (EOSAT/NASA, 
1986)  on  Sea-viewing,  Wide-Field-of-View  Sensor  (SeaWiFS). 

A  brief  description  of  what  is  known  about  the  atmospheric  corrections  and  sensor  degradation  is  given 
in  section  III.A  and  III.B  respectively. 

A.  Atmospheric  corrections 

The  majority  of  the  images  analyzed  in  the  initial  papers  on  the  CZCS  were  atmospherically  corrected 
using  Gordon’s  single  scattering  model  (Gordon  and  Clark,  1979).  Since  then  investigators  in  general,  have 
corrected  CZCS  data  using  an  improved  atmospheric  correction  algorithm  which  includes  multiple  scattering 
(Gordon  and  Castano,  1987;  Gordon  et  al.,  1988;  Gordon,  1990).  In  general  it  was  found  that  the 
atmospheric  contribution  to  the  upwelling  signal  detected  at  the  satellite  sensor  amounted  to  between  80- 
90%  of  the  total  signal. 

Thomas  and  Strub  (1990)  and  Strub  et  al.  (1990)  discussed  in  great  detail  a  West  Coast  Time  Series 
(WCTS)  of  the  California  Current  System  (CCS)  made  using  CZCS  data  from  mid  1979  to  mid  1986.  Due 
to  sensor  problems  that  occurred  in  1984  and  1985,  the  WCTS  did  not  include  data  for  these  years.  The 
atmospheric  correction  algorithm  used  to  process  the  WCTS  data  (single  scattering)  is  known  to 
overestimate  pigment  concentrations  at  large  solar  zenith  angles  present  during  winter  at  higher  latitudes 
(32°  N)  -  often  >3.0  mg  m3  or  greater  than  3  times  higher  than  found  during  the  1984  and  1985  winter 
cruises.  For  this  reason  the  months  of  November-February  of  each  year  were  eliminated  from  this  study.  As 
in  several  other  investigations,  the  estimate  of  the  total  error  in  comparing  CZCS  derived  pigment 
concentrations  with  the  in  situ  data  is  given  as  a  factor  of  2.  To  date  the  CCS  data  set  has  not  been 
analyzed  using  Gordon's  new  multiple  scattering  atmospheric  correction  model. 

In  a  recent  paper  Gordon  (1 990)  included  his  multiple  scattering  atmospheric  correcting  algorithm  along 
with  suggestions  on  changes  that  can  be  applied  to  this  algorithm  to  improve  the  accuracy  of  derived  optical 
parameters  from  future  ocean  color  sensors. 
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B.  Sensor  degradation 

Gordon  et  al.  (1983b)  showed  that  there  had  been  a  general  degradation  of  CZCS  sensors  starting  in 
1979.  They  were  able  to  quantify  this  degradation  by:  (1)  computing  the  water-leaving  radiance  for  imagery 
acquired  in  regions  where  the  water-leaving  radiance  was  known,  or  could  be  independently  estimated  and 
(2)  by  adjusting  the  sensor  calibration  to  force  agreement  between  the  two  radiances.  The  decay  of  the  blue 
band  at  443  nm  was  found  to  have  decreased  to  approximately  80%  of  its  initial  value  by  orbit  number 
20,000  (Sept  1982).  In  addition,  there  were  fairly  large  errors  associated  with  this  decay  factor  of  ±30-50%. 
Hovis  et  al.  (1985)  performed  aircraft  measurements  to  determine  the  degradation  of  the  CZCS  channels. 
They  were  able  to  confirm  the  major  source  of  degradation  was  due  to  loss  of  reflectance  from  the  large 
optics  that  were  exposed  to  the  spacecraft  environment,  namely,  the  scan  mirror  and  the  primary  and 
secondary  mirrors  of  the  telescope.  They  showed  that  the  response  of  the  blue  channel  (443  nm)  degraded 
approximately  25%  over  a  period  of  4  years  and  7  months  from  launch  to  May  1983.  The  degradation  in  the 
520  nm  band  was  3%. 

Mueller  and  Lange  (1989)  published  research  on  bio-optical  provinces  stating  that  there  existed  a  bias 
of  1 .5-2.0  between  the  CZCS  K490  data  and  in  situ  measurements.  The  primary  source  of  this  bias  error  was 
suspected  to  be  an  inaccurate  history  of  the  sensor’s  radiometric  sensitivity.  Various  investigators  have 
stated  their  interest  in  performing  a  full  7-year  CZCS  time  series  of  the  radiometric  consistency.  There  is 
evidence  that  not  only  did  the  sensitivity  of  the  sensors  degrade  over  the  years,  but  also  the  CZCS  changed 
its  sensitivity  at  an  irregular  rate  through  its  operating  lifetime. 

IV.  Preliminary  Satellite  Error  Analysis  Model 

The  analysis  of  CZCS  data  to  derive  phytoplankton  concentration  and  other  associated  variables  has 
been  discussed  in  sections  l-lll.  A  number  of  articles  have  pointed  out  that  the  atmosphere,  instrumental 
factors,  unresolved  (subpixel-sized)  clouds,  sensor  calibration  changes,  as  well  as  the  variation  of 
phytoplankton  concentrations  with  depth  below  the  ocean  surface,  all  contribute  to  the  total  error  budget 
of  satellite  derived  ocean  color  parameters.  In  addition,  errors  will  arise  when  comparing  data  acquired  on 
differ?^  dat"'-  (different  sense's,  or  altered  calibrations,  different  sun  angles,  different  atmospheres)  or 
different  scan  angles  of  the  same  or  different  images.  These  errors  will  be  bandpass  dependent,  which  result 
in  even  larger  errors  when  more  than  one  band  is  being  used  in  the  analysis.  For  instance,  using  a  two-band 
ratio,  as  used  in  the  present  analysis  to  derive  the  phytoplankton  concentration  (C)  or  the  diffuse  attenuation 
coefficient  (K),  the  error  will  be  approximately  double  the  error  in  one  band. 

The  generic  model,  which  is  given  in  this  section,  was  written  for  the  MS  -DOS  bar-ed  personal  computer. 
However,  more  extensive  programs  may  be  written  using  extended  memory  that  is  available  in  the  UNIX 
version  of  the  mathematical  software  used.  The  advantage  of  a  computer  model  is  that  many  potential  cases 
can  be  investigated  rapidly,  so  that  a  parametric  envelope  for  the  acquisition  of  useful  data  may  be 
developed,  resulting  in  a  realistic  appraisal  of  the  accuracy  of  the  derived  ocean  color  (phytoplankton 
concentration)  parameter.  In  this  way,  an  understanding  of  the  anticipated  precision  of  data  derived  from 
the  CZCS,  or  from  future  satellite-borne  or  airborne  sensors  may  be  developed. 
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The  upwelling  radiance  Lw(e,0)  at  the  satellite  sensor  is  given  by  equation  (4),  while  Figure  2  is  a 
schematic  diagram  for  a  simplified  two-dimensional  model  used  to  describe  the  factors  used  in  this  equation. 


Figure  2.  Overall  picture  of  interactions  between  incident  and  reflected  radiance,  the  atmosphere,  and  the 
Earth  s  surface. 


Input  Parameters  and  Definitions 

Calculations  of  upwelling  (water  leaving)  digital  radiance  LW($)  at  the  satellite  for  one  channel  is  given 
by  the  following  expression 

LW(4>)  =  sensor  gain  [solar  irradiance  X  atmospheric  transmission 

for  incident  and  reflected  paths  X  ocean  reflectance  +  atmospheric 

backscatter]  +  sensor  offset  (4) 


which  becomes 


LW(0)  =  GA.  (E(0)  .ri(9)  .rr(0)  .p  (0)  +  T(0))  +  OF,  (5) 


where  the  various  expressions  and  values  for  the  parameters  follow. 
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The  reflectance  of  seawater,  p  (<f>)  is  given  by  equation  (6),  where  a  4%  reflection  coefficient  at  500  nm 
has  been  assumed  for  seawater.  In  addition,  a  weak  view  angle  dependence  has  been  included. 


P  (n) 


0.04-  exp 


tt  y 

0.1  cos 

<*>(  n) - 

18  0J 

(6) 


note:  4>(n)  =  5(n-l2)  -  view  zenith  angle  in  degrees,  and  degrees  are  converted  to  radians 
n  =  0...24  Counter  (running  variable) 

The  atmospheric  transmission  along  the  path  from  the  Sun  to  the  Earth's  surface  is  given  by  equation 

(?) 


t i(0)  =  exp 


-a 

n 

cos 

6  - 

180 

. 

(7) 


where  ©  is  the  solar  zenith  angle  In  degrees,  and  a  Is  the  extinction  coefficient,  which  has  been  assumed 
to  be  constant  over  the  path  length  and  equal  to  0.35. 


The  atmospheric  transmission  along  the  path  from  the  Earth’s  surface  to  the  sensor  is  given  by  equation 

(8) 


tr(n)  =  exp 


-a 

n 

cos 

<P(  n) - 

180 

_ 

(8) 


The  approximation  for  the  backscatter  r  (n)  at  the  sun  zenith  angle  (0)  and  view  angle  (4>)  is  given  by 
equation  (9). 


" 

ir 

T(n)  =  E (©)  p  (n)  exp 

-b 

1  +  cos 

I©  -  0(n) ! - 

180 

0) 

B 


where 


E(0) 


Isol- cos 


e 


180 


Irradiance  on  sea 


(10) 


Isol  =  Solar  irradiance  above  atmosphere  for  0.1  micron  wide  bandpass  dose  to  500  nm: 
expressed  in  W/m2  (was  set  equal  to  135) 

b  =  constant,  which  is  set  to  yield  the  approximate  amount  of  backscatter.  It  was  set  equal  to  0  4  in 
the  present  calculations  (amounting  to  approximately  50%  backscatter). 

OF  =  sensor  offset  (=  2) 

GA  =  Sensor  gain  for  8  bit  quantization  (=  80) 


Consider  the  case  where  the  upwelling  (water  leaving)  radiance  at  nadir  view  angle  (4>  =  0  (n  =  12))  and 
solar  elevation  (0)  of  -45.  Equation  (5)  becomes 


LW(n)  =  GA  (E(-45)  ti (-45)  rr(12)  p  (12)  +T(12))  +  '  *F  (11) 


Assume  now  a  change  of  3%  error  (el)  in  the  sensor  offset  (OF)  and  a  10%  error  e2  in  the  sensor  gain 
(GA).  Equation  (11)  now  is  changed  to  read 


LWe(n)  =  GA 


1 


(E(-45) • r i ( —4 5 )  rr(n)  p  (n) 


+  T(n)) 


+ 


1 


+ 


e2 


100 


■  OF 


(12) 


The  percent  change  in  upwelling  radiance  with  introduction  of  sensor  errors  and  sun  angle  and  scan 
angle  effects  for  one  channel  is  given  by  equation  (13) 


ER(n) 


100- 


LWe(n)  -  LW(n) 
LW(n) 


(13) 
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The  total  error  in  rationing  two  channels  is  approximately  double  that  in  one  channel. 

The  following  calculations  shown  in  Table  1  and  Figure  3  give  the  difference  from  the  radiance  from 
nadir  for  prelaunch  calibration  (gain  and  offset)  due  to  varying  view  zenith  angle  and  to  errors  in  sensor  gain 
and  offset.  Here  a  solar  zenith  angle  of  45°  is  assumed.  Change  in  solar  zenith  angle  and  in  atmospheric 
transmission  may  yield  even  larger  errors. 

These  calculations  are  given  only  as  representative  of  the  type  of  information  that  can  be  derived 
relatively  quickly.  Such  a  modeling  technique  can  be  used  to  perform  an  absolute  error  analysis  of  each 
source  of  satellite  sensor  error.  The  result  is  a  total  error  budget  for  the  system. 


Table  1 .  Calculations  of  radiances  and  errors  associated  for  various  view  zenith  angles. 


Recorded  radiance  difference 

due  to  view  zenith  angle  and  Resulting  error  in 
View  zenith  angle  to  sensor  calibration  errors  radiance  ratio 


0(n)  LWe(n)  -  LW(n)  2  ER(n) 


-60 

-47.675 

-30.025 

-55 

-36.796 

-23.174 

-50 

-27.539 

-17.343 

-45 

-19.452 

-12.25 

-40 

-12.221 

-7.696 

-35 

-5.624 

-3.542 

-30 

0.497 

0.313 

-25 

6.253 

3.938 

-20 

11.726 

7.385 

-15 

16.97 

10.688 

-10 

22.022 

13.869 

-5 

26.902 

16.943 

0 

31.617 

19.912 

5 

36.159 

22.772 

10 

40.509 

25.512 

15 

44.633 

28.11 

20 

48.482 

30.533 

25 

51.986 

32.74 

30 

55.055 

34.673 

35 

57.569 

36.256 

40 

59.368 

37.389 

45 

60.246 

37.942 

50 

59.924 

37.74 

55 

58.027 

36.545 

60 

54.039 

34.033 
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(a)  Modeled  water  reflectance  as  a 
function  of  view  zenith  angle 
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(b)  Modeled  atmospheric  transmission 
as  a  function  of  view  zenith  angle 
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(c)  Modeled  backscattered  radiance 

as  a  function  of  view  zenith  angle 
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(d)  Modeled  recorded  digital  radiance 
from  a  sensor  with  calibration 
error,  shown  as  afunction  of  view 
zenith  angle 
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Figure  3.  Results  of  calculations  of  various  parameters  affecting  the  final  derived  radiance  at  the  satellite. 
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Image  based  studies  will  lead  to  a  knowledge  of  the  variance  of  recorded  radiances  over  regions  in  each 
bandpass.  Appropriate  theoretical  studies  will  advance  our  understanding  of  the  impact  of  such  variations 
on  the  derived  pigment  concentrations  (C).  Further,  image  band  and  modeling  studies  will  result  in  an 
improved  understanding  of  the  impact,  for  instance,  of  the  atmosphere  on  the  derived  values  for  C.  The 
result  is  a  more  accurate  assessment  of  the  accuracy  of  the  derived  value  for  C  under  a  variety  of 
bandpasses  and  imaging  conditions.  In  this  way  imaging  and  analysis  may  be  optimized  to  most  accurately  ' 
map  C.  In  addition,  studies  investigating  the  impact  of  different  imaging  averaging  procedures  need  to  be 
performed,  since  quite  often  the  procedure  currently  followed  is  to  use  one  pixel  in  25  as  representative  of 
the  radiance  instead  of  averaging  the  total  of  25  pixels  in  the  region.  It  is  well  known  that  the  accuracy  of 
index  maps,  using  LAND  SAT  satellite  data,  is  greatly  impacted  by  the  averaging  procedures  employed. 

V.  Conclusions  and  Recommendations 

In  summary,  there  are  a  large  number  of  sources  of  error  that  one  must  consider  when  using  the  C2CS 
images  to  derive  optical  properties  of  the  ocean.  On  a  global  basis  there  is  general  agreement  that  the 
pigment  concentration,  which  is  directly  related  to  the  diffuse  attenuation  coefficient  (K)  can  be  determined 
from  CZCS  data  only  to  within  a  factor  of  2,  in  comparison  with  in  situ  measurements.  In  a  much  more 
regional  or  restrictive  area,  there  is  agreement  that  this  error  may  be  reduced  to  approximately  +30-40%. 

In  order  to  reduce  the  error  to  this  magnitude  extreme  care  on  making  appropriate  corrections  must  be 
made  on  a  pixel -to-pixel  basis.  Only  a  few  reports  give  errors  that  are  this  small. 

Of  particular  importance  to  this  summary  report  are  the  papers  by  Strub  et  al.  (1990)  and  Thomas  and 
Strub  (1990),  which  were  mentioned  in  sections  II  and  III  to  include  a  complete  time  series  of  the  CZCS  for 
the  California  Current  data.  A  close  look  at  this  reference  and  the  techniques  used  in  this  analysis  should 
be  of  real  value  to  additional  time  series  analysis  performed  on  the  CZCS  data  at  the  Naval  Oceanographic 
and  Atmospheric  Research  Laboratory  (NOARL). 

It  is  our  recommendation  that  a  modeling  program  be  initiated  at  NOARL  to  investigate  the  total  error 
budget  associated  with  using  the  CZCS  data  to  derive  parameters  relating  to  the  surface  optical  properties 
of  the  ocean  waters  for  the  various  ADI  areas.  In  section  IV  a  generic  model  for  a  PC  was  discussed  which 
could  be  used  to  analyze  the  total  error  budget.  It  is  emphasized  that  this  is  strictly  a  generic  model,  which 
by  necessity  will  become  much  more  complicated. 
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